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MEASUREMENT  OF  HYDROACOUSTIC  PARTICLE  MOTION  BY  HOT-FILM  ANEMOMETRY 


INTRODUCTION 

Most  hydroacoustic  detectors  are  sensitive  to  pressure  fluctuations. 

When  the  acoustic  field  Is  simple  and  of  known  signature,  one  may  compute 
particle  velocity  from  a  single  local  pressure  measurement.  For  more  compli¬ 
cated  fields  it  is  necessary  to  map  the  pressure  field  in  both  magnitude  and 
phase  in  order  to  derive  the  particle  acceleration  by  spatial  differentiation. 
It  is  attractive  to  measure  the  particle  motion  independently,  as  a  check  on 
the  computations  and  as  an  alternate  method  of  acoustic  measurement.  Examples 
of  such  measurements  are  given  in  References  1  and  2. 

The  method  of  hot-film  or  hot-wire  anemometry  offers  a  new  approach  to 
the  detection  of  particle  motion  in  acoustic  fields  in  liquids  and  gases  [3]. 
It  determines  the  amount  of  heat  transfer  from  an  electrically  heated  film  or 
wire  to  the  surrounding  fluid  medium.  This  heat  transfer  depends  on  a  number 
of  factors,  but  under  the  proper  conditions  it  is  predominantly  a  function  of 
the  kinematic  parameters  of  the  fluid  particles  in  the  medium.  Review 
articles  on  this  technique  are  found  in  References  4,  3,  and  6. 

The  method  of  hot-film  anemometry  has  found  application  wherever  its  fast 
response  time  and  small  physical  size  are  of  advantage.  Thus  it  has  become 
the  measurement  device  of  choice  in  the  study  of  turbulence.  To  the  best  of 
the  author's  knowledge,  no  other  reports  on  application  of  hot-film  anemometry 
to  hydroacoustics  are  found  in  the  literature.  A  computer  search  by  the 
Defense  Technical  Information  Center  did  not  produce  any  relevant  articles. 
There  is  extensive  literature  on  the  influence  of  intense  acoustic  fields  or 
of  the  vibration  of  a  given  body  on  its  heat  transfer  characteristics 
[7-12],  but  no  mention  is  made  of  the  principle  of  using  the  heat  transfer  to 
measure  the  acoustic  radiation  for  its  own  sake. 
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It  appears  to  be  tacitly  assumed  in  practically  all  applications  of  hot- 
film  anemometry  that  the  operation  of  the  device  may  be  described  by  the 
precepts  of  the  forced  convection  regime.  As  a  consequence  the  calibration 
relations  are  cast  in  the  form  of  a  relation  between  output  voltage  and 
imposed  velocity,  usually  as  a  variation  on  King's  law  [13].  Two  aspects  of 
acoustic  particle  motion  militate  against  application  of  the  same  type  of 
relation  in  this  case.  First,  the  particle  speeds  are  quite  small  for 
acoustic  fields  of  modest  intensity  so  that  free  convection  is  not  negligible. 
Secondly,  the  periodicity  implies  that  the  same  fluid  will  repeatedly  pass  by 
the  sensor  with  only  incomplete  cooling  by  the  ambient  fluid  during  one 
period. 

A  theoretical  analysis  of  the  flow  and  heat  transfer  phenomena  in  the 
presence  of  an  acoustic  field  appears  quite  formidable.  In  order  to  create  a 
scheme  for  arranging  calibration  measurements,  use  is  made  of  the 
dimensionless  numbers  that  appear  when  a  dimensional  analysis  is  made  of  the 
equations  of  motion  and  heat  diffusion. 

The  existence  and  importance  of  free  convection  lead  to  various  modes  of 
operation  of  hot-film  anemometry  in  hydroacoustics.  When  the  acoustic  veloc¬ 
ity  is  parallel  to  gravity,  the  free  convection  acts  as  a  dc  bias  flow  for  the 
acoustic  velocity.  Therefore,  the  anemometer  response  is  expected  and  found 
to  have  a  frequency  equal  to  the  acoustic  frequency  and  an  amplitude  propor¬ 
tional  to  the  particle  motion.  When  the  acoustic  velocity  is  horizontal,  it 
is  not  accompanied  by  a  dc  bias  flow.  Because  of  the  symmetry  of  the  physical 
situation  for  the  movement  of  the  particles  to  one  side  or  the  other  of  the 
sensor ,  one  expects  an  anemometer  response  with  twice  the  frequency  of  the 
acoustic  field  and  an  amplitude  proportional  to  the  square  of  the  particle 
motion. 

A  third  mode  of  operation  of  the  anemometer  may  be  induced  by  creating  an 
imposed  dc  bias  flow  by  a  nozzle  that  directs  a  jet  of  fluid  onto  the 
sensor.  The  sensitivity  of  the  instrument  is  thereby  increased,  and  one  is 
largely  free  from  the  dependence  on  the  direction  of  gravity.  The  direction 
of  this  jet  determines  a  preferred  direction  which  improves  the  directional 
properties  of  the  anemometer.  A  second  directional  sensitivity  is  supplied  by 
the  axis  of  the  cylindrical  sensor:  the  preferred  direction  for  the  flow  is 
perpendicular  to  this  axis  for  all  three  modes  of  operation. 
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The  next  section  offers  a  theoretical  analysis  and  a  formulation  of  the 
pertinent  dimensionless  numbers.  The  creation  and  calibration  of  flow  fields 
needed  to  test  the  anemometer  are  described,  and  the  data  are  given  as  a 
function  of  the  relevant  dimensional  parameters.  A  comprehensive  view  of  the 
data  is  given  by  empirical  relations  in  terms  of  dimensionless  numbers.  The 
computation  of  the  various  dimensionless  numbers  is  discussed  in  Appendix  A. 
Appendix  B  describes  the  construction  of  a  trough  and  shaker  arrangement  used 
to  produce  horizontal  particle  motion. 

1.  THEORETICAL  ANALYSIS 


1.1.  Basic  Equations  and  Processes 

The  fundamental  equations  that  govern  the  phenomenon  of  heat  transfer  in 
hot-film  anemometry  are  the  conservation  of  mass,  the  equation  of  motion,  and 
the  conservation  of  energy.  Usually  these  equations  are  presented  in  a 
simplified  form,  known  as  the  Boussinesq  approximation  [14,15].  The 
continuity  equation 

|f  +  $  *  (Pu)  -  0,  (1) 

(where  p  is  the  density  of  the  medium  and  u  is  the  particle  velocity  vector) 
is  simplified  to  the  incompressibility  relation 

$  •  u  »  0.  (2) 


The  suppression  of  the  time  derivative  of  the  density  even  in  the  presence  of 
acoustic  radiation  is  based  on  the  fact  that  the  size  of  the  sensor  is  small 
compared  with  the  wavelength  of  sound  in  the  region  studied. 

The  equation  of  motion  in  this  approximation  is 

|f  +  (u  •  $)  u  ■  -  ~  $p  +  gaCit  +  vV2u ,  (3) 

o 

where  a  -  (1/V)(3V/3T)  is  the  thermal  expansion  coefficient  and  V  is  volume; 

P 

0  ■  T  -  T  ,  the  difference  between  the  local  temperature  and  the 
temperature  T0  of  the  undisturbed  medium; 


v  «*  kinematic  viscosity; 

pQ  *  density  of  the  undisturbed  medium  (constant); 
g  *  acceleration  of  gravity; 
ic  “a  unit  vector  vertically  upward;  and 
p  *  dynamic  pressure. 

The  conservation  of  energy  reduces  to  a  temperature  equation 

+  (u  •  $)T  -  <72T,  (4) 

where  <  is  the  thermal  diffusivity  equal  to  k/pcp,  with  k  the  heat  conduction 

coefficient  and  cp  the  specific  heat  at  constant  temperature.  The 

approximation  implied  by  Eqs.  (3)  and  (4)  amounts  to  setting  all  physical 

parameters  of  the  medium  constant  except  the  density  insofar  as  it  causes 

buoyancy.  The  density  difference  in  the  buoyancy  term  is  approximated  by  a 

linear  expression  in  the  temperature  or  2£.  =  -a9  • 

o 

The  measurement  of  acoustic  particle  motion  by  means  of  hot-film 
anemometry  consists  of  two  interacting  processes.  The  heat  transfer  from  the 
hot  film  in  a  quiescent  fluid,  in  addition  to  the  ever-present  heat 
conduction,  is  due  to  convection  by  the  buoyancy  of  the  fluid  that  is  heated 
by  the  film  itself.  This  process  is  called  free  convection.  The  second 
process  is  that  of  an  incompressible  fluid  motion  about  the  cylindrical  sensor 
accompanying  the  acoustic  radiation. 

The  first  process  may  be  described  by  Eqs.  (2),  (3),  and  (4),  where  the 
time  derivatives  3u/3t  and  3T/3t  are  omitted,  since  the  fluid  motion  near  the 
cylinder  is  steady  and  becomes  unstable  only  several  diameters  away  from  the 
cylinder.  The  equation  of  motion,  Eq.  (3),  and  the  temperature  equation, 

Eq.  (4),  are  coupled  and  have  to  be  solved  simultaneously.  This  cannot  be 
done  analytically,  even  for  simple  geometries;  one  has  to  rely  on  empirical 
relationships  [16-18].  The  second  process,  flow  past  a  cylinder  perpendicular 
to  the  flow,  is  more  amenable  to  analytical  or  numerical  treatment. 
Unfortunately,  even  with  more-or-less  extensive  knowledge  of  the  two  separate 
processes  it  is,  in  general,  not  possible  to  add  their  effects  because  of  the 
nonlinear  advective  acceleration  (u  •  ^)u.  In  case  this  term  is  small 
compared  with  local  acceleration  3u/3t,  one  may  add  the  two  solutions  for  the 
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free  convection  and  imposed  acoustic  flow.  One  might  intuitively  feel  that 
the  horizontal  particle  motion  due  to  the  acoustic  field  is  reasonably 
independent  of  the  motion  due  to  the  thermal  plume  in  the  free  convection  and 
could  be  independently  solved,  as  if  it  were  forced  convection.  The  resulting 
velocity  field  could  then  be  introduced  into  the  temperature  equation  and 
solved  for  the  relevant  heat  transfer. 

1.2.  Dimensional  Analysis 

The  complexity  of  the  equations  governing  the  1  *ransfer  by  a  hot 
film,  in  the  presence  of  acoustic  radiation,  renders  ’  analytical  or 
numerical  solution  a  formidable  task.  In  order  to  establish  a  calibration  of 
the  hot-film  anemometry  to  serve  as  a  tool  in  acoustic  research,  one  may 
utilize  the  concepts  of  dimensional  analysis  to  reduce  the  number  of  variables 
influencing  the  phenomenon  to  a  smaller  set  of  dimensionless  groups.  For 
reference,  Eqs.  (3)  and  (4)  are  repeated  here,  with  a  short  description  of  the 
meaning  of  each  term  written  under  the  equation. 


Equation  of  Motion: 


3u 

TT 

+  (u  •  hn  = 

♦ 

o 

ga0lt  + 

w2lr  (5) 

local 

advective 

pressure 

buoyancy 

viscosity 

acceleration 

acceleration 

force 

Temperature  Equation: 

3T 

W  + 

(it  •  $)T 

*  icV2T 

local 

advection 

conduction 

heating  rate 

of  heat 

of  heat 

Dimensional  analysis  consists  in  comparing  the  order  of  magnitude  of  the 
various  terms.  The  ensuing  dimensionless  groups  serve  as  the  independent 
variables  in  the  sought-for  empirical  relations. 
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The  relevant  dependent  variable  is  the  heat  flux  H,  expressed  in 
dimensionless  form  by  the  Nusselt  number  Nu  =  (Hd/kOg),  where  d  is  the 
diameter  of  the  sensor  and  0g  =  Tg  -  TQ  is  the  difference  in  temperature 
between  sensor  and  undisturbed  medium. 

1.3.  Forced  Convection  Versus  Free  Convection 

An  important  distinction  in  heat  transfer  by  fluids  is  that  between 
forced  convection  and  free  convection.  In  the  first  case  there  is  an  imposed 
flow  of  sufficient  strength  to  make  the  buoyancy  term  negligible.  The 
equation  of  motion  then  is  not  coupled  to  the  temperature  equation  and  the 
equation  of  motion  may  be  independently  solved.  The  solution  for  the  velocity 
field  is  entered  into  the  temperature  equation  through  the  advected  heat 
term  u  •  VT,  and  this  equation  is  then  solved  for  the  temperature.  In  free 
convection  the  two  equations  have  to  be  solved  simultaneously,  since  the  term 
driving  the  motion  is  the  buoyancy  iti ’If.  For  the  present  application  the 
distinction  is  important  since  one  likes  to  know  what  the  value  of  an  imposed 
bias  flow  may  be  in  order  to  qualify  for  the  designation  of  forced  convection. 

Suppose  that  a  bias  flow  with  typical  speed  U  is  imposed  on  the  sensor. 

To  judge  its  significance  one  compares  the  buoyancy  term  ga0g  in  Eq.  (5)  with 
either  the  advective  acceleration  term  (u  •  $)u  or  the  viscosity  term  vV2u  , 
depending  on  which  of  these  is  dominant.  The  relevant  dimensionless  group  to 
compare  advective  acceleration  (inertia)  with  viscosity  is  the  Reynolds  number 
Re  =  Ud/v.  If  Re>>l ,  one  compares  buoyancy  with  inertia,  which  results  in  the 
Richardson  number  Ri  =  ga0d/U2.  If  Re«l ,  one  compares  buoyancy  with 
viscosity  and  the  pertinent  number  is  ga0d2/Uv,  equal  to  RiRe  (no  separate 
name  was  found  for  this  dimensionless  group).  Thus  forced  convection  is 
present  if  Re»l  and  Ri«l ,  or  if  Re<<l  and  ReRi<<l.  Free  convection  prevails 
if  Re»l  and  Ri»l ,  or  Re<<l  and  ReRi>>l.  Any  intermediate  case  is  called 
mixed  convection.  The  borderline  speed  Ucf  is  obtained  by  setting  either 
Ri  *  1  and  solving  for  U,  or  by  setting  ReRi  =  1,  depending  on  whether  this 
given  Ucr  makes  Re>l  or  Re<l ,  respectively. 

In  Table  1  the  values  of  Ucr  are  given  for  the  three  sensor  diameters  and 
three  temperature  differences  between  sensor  and  medium  used  in  this  study. 
Both  water  and  ethylene  glycol  were  studied.  Except  for  the  point  marked  with 
an  asterisk,  viscosity  dominates  over  inertia.  In  Table  2  values  for  U  are 


given  that  make  the  Reynolds  number  equal  to  1.  The  importance  of  this  table 
is  that  one  may  judge  for  which  values  of  1)  the  convection  is  of  the  forced 
type,  while  still  Re<l.  The  latter  condition  defines  a  flow  regime  known  as 
creeping  flow.  For  such  flow  about  a  cylinder,  Stokes  has  given  an  analytic 
solution,  refined  by  Oseen  [19,20].  Thus,  in  some  cases  of  an  imposed  flow 
with  the  proper  speed,  one  may  use  this  analytic  solution,  insert  it  into  the 
temperature  equation,  and  find  a  solution  of  the  latter,  presumably  in  numeric 
form.  Usually,  though,  one  has  to  settle  for  an  empirical  relation  for  the 
given  heat  transfer  in  terms  of  the  relevant  dimensionless  number(s). 

For  heat  transfer  from  a  cylinder,  this  is  usually  designated  as  "King's  law”, 
after  the  author  of  the  first  investigation  into  this  subject  [13]. 

Morgan  [21]  proposes  a  relation  of  the  form 

Nu  =  (C  +  D  Re11)  PrP ,  (7) 

where  Pr  is  the  Prandtl  number  Pr  =  v/k,  and  lists  values  for  the  constants  C, 
D,  n,  and  p  from  many  sources.  Blackwelder  [22]  writes  the  relation  in  a  more 
general  form 

Nu  =  [A(Pr,aT)  +  B(Pr,aT)Ren]  (1  +  aT/2)m,  (8) 

where  a^.  is  the  overheat  ratio  aT  =  0g/To  (TQ  in  Kelvin). 

The  aspect  ratio  2/d,  where  l  is  the  length  of  the  sensor,  would  also  be 
expected  to  be  important.  This  is  especially  true  for  hot  films,  where  this 
ratio  is  much  smaller  than  for  wires,  and  thus  the  conditions  are  farther 
removed  from  the  ideal  case  of  an  infinite  cylinder. 

This  proliferation  of  dimensionless  numbers  and  relations  was  addressed 
by  Blackwelder  [22]  in  a  thoughtful  comment: 

"It  is  distressing  (especially  to  the  uninitiated)  that 
there  is  such  a  lack  of  agreement  concerning  the  basic 
equation  of  hot-wire  anemometry. . .  These  variations  and 
deviation  from  an  exact  universal  equation  necessitate 
that  each  hot-wire  and  hot-film  be  calibrated 
individually. .. before  it  can  be  used  to  record  and 
interpret  data." 
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This  quote  may  warn  against  too  optimistic  a  hope  of  finding  a  comprehensive 
empirical  equation  for  axl  circumstances  of  sensor  type,  frequency, 
temperature,  and  medium  in  the  application  to  hydroacoustics. 

The  Eqs.  (7)  and  (8)  may  be  used  to  relate  particle  speed  to  measured  dc 
voltage  in  the  case  of  imposed  bias  flow. 


Table  1  -  Critical  speed  Ucr  in  mm/s  as  a  function  of  sensor 

diameter  d  and  temperature  difference  ©g.  For  U»Ucr 
forced  convection  obtains,  for  U«Ucr  free  convection. 
The  values  for  Ucr  are  found  from  ReRi  =  1,  except  the 
point  marked  *,  where  Ri  *  1. 
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- > 

Water 
d  in  pm 

Ethylene  glycol 

25 

50 

150 

25 

50 

150 

0.039 

0.15 

mm 

0.0051 

0.021 

0.19 

0.11 

0.43 

3.9 

0.014 

0.057 

0.51 

0.21 

0.85 

5.4* 

0.029 

0.11 

1.0 

Table  2  -  Speed  U  in  mm/s  for  which  Re  =  1. 

For  Re<l  one  may  apply  the  results  of 
Stokes  and  Oseen  for  creeping  flow  [19,20] 


Water 

— >  d  in  pm 

Ethylene  glycol 

25 

50 

150 

25 

50 

150 

30.8 

15.4 

5.13 

474 

237 

79 

26.3 

13.2 

4.39 

348 

174 

58 

22.8 

11.4 

3.81 

258 

129 

A3 

1.4.  Free  Convection 

By  definition  the  term  free  convection  implies  that  the  buoyancy 
term  go©  dominates  the  equation  of  motion.  The  Reynolds  number  still  would  be 
important  to  distinguish  between  inertia-dominated  and  viscosity-dominated 
flow  regimes.  Obviously  Re  cannot  be  directly  evaluated  since  the  value  of  U 
is  not  known  a  priori.  It  has,  in  turn,  to  be  estimated  from  dimensional 
considerations.  The  buoyancy  term  must  be  matched  by  either  advection  or 
viscosity  in  the  equation.  If  buoyancy  and  advection  match,  the  value  of  U  is 
U  *  (ga0d)^2  and  inserting  this  U  into  the  Reynolds  number  gives  a  group 
[ (ga0d3)/v2] ^  •  If  buoyancy  and  viscosity  match,  one  has  U  «  (ga0d2)/v  and 
the  Reynolds  number  is  transformed  into  (ga©d3)/v2.  The  Grashof  number  is 
defined  by  Gr  =*  (ga©d3)/v2  and  one  sees  that  it  takes  the  place  of  the 
Reynolds  number:  inertia  dominates  for  large  Grashof  number,  viscosity 
dominates  for  small  Gr,  with  the  understanding  that  the  situation  is  one  of 
free  convection  where  the  buoyancy  term  is  the  prime  mover. 

A 

One  may  apply  the  same  reasoning  to  the  Peclet  number,  which  compares 

heat  advection  with  heat  conduction.  It  follows  that  the  product  GrPr  takes 

the  place  of  the  Peclet  number  for  small  Gr;  for  large  Gr  it  is  the 

combination  Gr  ^2  Pr.  The  product  GrPr  is  also  known  as  the  Rayleigh  number, 

3 

Ra  =  GrPr  =  (gaGd  )/vtc. 

An  empirical  relation  connecting  heat  transfer  to  the  independent 
variables  in  free  convection  is  given  by  Morgan  [23] ,  who  lists  values  for  the 
constants  A,  B,  m  from  various  sources,  as 

Nu  =  A  +  B  (GrPr)m.  (9) 

This  relation  is  used  to  interpret  dc  heat  transfer  in  Section  3.1.  Jaluria 
[24]  gives  a  relation  Nu  =  C(GrPr)^  where  C  is  a  function  of  the  Prandtl 
number. 

1.5.  Response  to  Acoustic  Fields 

It  was  assumed  in  the  foregoing  sections  that  the  phenomena  were  steady. 
The  Instabilities  in  the  thermal  flow  were  ignored,  since  they  appear  only 
several  diameters  away  from  the  sensor.  If  one  introduces  a  harmonic  acoustic 
particle  motion  with  angular  frequency  u>,  the  time  derivatives  in  Eqs.  (5)  and 


(6)  are  to  be  reckoned  with.  The  first  dimensional  estimate  concerns  the 
ratio  of  advective  acceleration  to  the  local  derivative,  resulting  in  the 
dimensionless  number  U/md,  where  U  is  the  rras  value  of  a  harmonic  vibration. 

If  one  introduces  the  rms  displacement  vector  t  with  components  £,  n,  5,  this 
may  also  be  written  as  C/d  and  c/d  for  horizontal  and  vertical  acoustic 
motions,  respectively. 

For  horizontal  particle  motion  there  is  symmetry  for  the  motion  to  the 
right  or  to  the  left  of  the  film.  Therefore,  one  expects  that  the  frequency 
of  the  anemometer  response  is  twice  the  frequency  of  the  acoustic  field. 
Moreover,  the  first  term  in  an  expansion  of  the  response  in  powers  of  relative 
displacement  would  be  quadratic.  Thus,  for  modest  values  of  the  relative 
displacement,  one  expects  that  the  Nusselt  number  varies  harmonically  with 
twice  the  frequency  of  the  field,  and  that  the  rms  value  of  the  fluctuating 
component  NuflC  is  given  by 


Nu  =■  b(C/d)2 ,  (10) 

3C 

where  the  coefficient  b  is  independent  of  the  displacement  5.  Nu  is 
proportional  to  the  rms  value  of  the  ac  voltage  Vac-  For  the  vertical 
particle  motion  of  the  acoustic  field  there  is  a  dc  bias  flow  provided  by  the 
free  convection.  Therefore,  one  expects  the  frequency  of  the  response  to  be 
the  same  in  this  case  as  that  of  the  acoustic  field  and  the  Nusselt  number 
variation  to  be  linear  in  the  displacement  for  modest  values  of  this 
parameter.  Thus , 


Nu  -  a(c/d),  (11) 

ac 

where  the  coefficient  a  is  independent  of  the  displacement  £.  The  coef¬ 
ficients  a  and  b  might  be  functions  of  the  dimensionless  numbers  defining  the 
free  convection  field  Gr  and  Pr.  It  is  found  experimentally  that  the  coef¬ 
ficient  a  is  independent  of  frequency  within  the  range  covered  by  the  experi¬ 
ments;  the  coefficient  b,  though,  is  dependent  on  frequency.  Inspection  of 

the  temperature  equation  indicates  that  the  relevant  dimensionless  number  is 

o 

the  Fourier  number  Fo  ■  </u>d  .  It  compares  the  heat  carried  away  by  conduc¬ 
tion  with  the  local  variation  of  stored  heat  in  a  fluid  particle.  Another 

2 

frequency-dependent  dimensionless  number  is  the  Stokes  number  St  -  v/ud  , 
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which  compares  the  viscosity  with  the  local  acceleration.  Of  course,  one 
has  St  -  FoPr.  As  in  the  case  of  steady  flow,  keep  in  mind  the  aspect 
ratio  4/d  as  another  variable  able  to  influence  the  anemometer  response. 
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2.  EXPERIMENTAL  METHODS  AND  RESULTS 


The  experiments  conducted  with  the  hot-film  anemometer  fall  into  three 
categories.  The  response  to  horizontal  particle  motion  is  measured  in  a  newly 
developed  trough  and  shaker  arrangement  (see  Appendix  B).  For  vertical 
particle  motion,  the  NRL-USRD  type  G40  calibrator  is  used  [25].  The  velocity 
field  in  the  two  devices  is  calibrated  by  an  NRL-USRD  type  F61  standard 
hydrophone  [25].  The  pressure  field  is  determined  as  a  function  of  location; 
spatial  differentiation  then  gives  the  acceleration  according  to  the  relation 
(Euler's  equation) 


3u 

Tt 


(12) 


The  third  type  of  measurements  consists  of  the  response  to  acoustic  motion,  in 
the  presence  of  a  bias  flow  produced  by  a  nozzle,  in  the  G40  calibrator. 

Figure  1  depicts  the  experimental  arrangement.  The  hot-film  anemometer 
is  a  Thermo Systems  Inc.  System  No.  1050-1.  The  various  sensors  used  in  the 
experiment  are  shown  in  Fig.  2.  The  sensor  type  is  designated  by  a  number  and 
letter  combination;  e.g. ,  1210-20W.  Notice  that  the  last  two  digits  before 
the  dash  refer  to  the  shape  of  the  probe;  the  letter  W  indicates  a  quartz- 
coated  film  for  use  in  conducting  liquids.  The  two  digits  after  the  dash 
indicate  diameter:  10  corresponds  to  25  p,  20  to  50  p,  and  60  to  150  p.  The 
oscillator  is  an  HP3300A  function  generator  with  a  3305A  sweep  plug  in.  Its 
signal  drives  an  Optimation  Model  PA250M  power  amplifier  in  constant  current 
mode,  which  provides  the  power  for  driving  the  coil  of  the  GAO  calibrator  or 
the  shaker.  A  frequency  reference  signal  is  taken  from  a  standard  resistor  in 
series  with  the  coil  and  is  entered  into  the  reference  input  of  the  PAR  Model 
5204  Lock-In  Analyzer.  It  is  important  that  this  reference  signal  not  deviate 
too  much  from  the  optimum  value  of  1  V.  The  lock-in  analyzer  has  a  vector 
option  to  provide  magnitude  and  phase  of  the  anemometer  output,  the  phase  thus 
is  relative  to  the  phase  of  the  current  through  the  coil.  The  hot-film 
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anemometer  output  is  viewed  on  an  HP1200B  oscilloscope,  which  shall  be 
disconnected  when  measurements  are  taken  with  the  F61  hydrophone  in  order  not 
to  shunt  its  high  input  impedance.  The  lock-in  amplifier  has  a  If  -  2f 
setting  so  that  the  analyzed  frequency  may  be  taken  as  equal  to  the  field 
frequency  for  vertical  motion  and  twice  the  acoustic  frequency  for  horizontal 
motion.  The  output  signal  from  the  lock-in,  magnitude  or  phase,  is  plotted  on 
an  HP7015B  X-Y  recorder.  The  x  coordinate  is  derived  from  the  output  of  the 
sweep  generator  connected  with  the  oscillator. 


G40  OR 
SHAKER 


Fig.  1  -  Block  diagram  for  experimental  arrangement 


_ Configuration _ 

Model  1210  Standard  Straight  Probe 


Type 


-  VM?5l  Oa 


Size 


{•In 


d  =  50  V 


Model  1211  Sensor  Parallel  to  Probe 
Cross  Flow 


_  V  J 


Model  1212  Probe  with  90*  Bend 


5.  *  1  mm 

d  =  50  U 


£  ■  2  mm 


d  -150  u 


Model  1260 A  Straight  Probe-Miniature 


-  »n£0>  — M 

f— 1^8  1.625)  — H 


L  -  0.5 


1.5 

(.060)  01a 


Legend:  mm  (In.) 


3.2 

(.125)  01a 


d  *  25  P 


Figure  provided  by  courtesy 
of  TSI  Incorporated. 


Fig.  2  -  Sensors  used  in  experiments 
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2.1.  Horizontal  Particle  Motion 

A  trough  and  shaker  arrangement  was  especially  developed  to  create  a 
horizontal  periodic  particle  motion.  It  is  described  in  Appendix  B. 

Although  the  measurements  are  performed  in  the  center  of  the  trough, 
where  theoretically  no  vertical  motion  is  present,  a  1211  type  sensor  was  used 
in  the  earlier  measurements.  This  sensor  is  parallel  to  the  probe  and  probe 
holder,  and  therefore  is  parallel  to  the  direction  of  the  free  convection.  It 
is  expected  that  any  spurious  If  signals  would  be  attenuated  this  way,  by 
virtue  of  the  lesser  sensitivity  of  the  vertical  sensor  for  vertical  motion 
(see  the  discussion  of  vertical  particle  motion),  by  about  a  factor  of  four. 

It  was  found  that  no  2f  signal  is  produced  by  whatever  spurious  signals  are 
present.  Thus  the  sensors  with  diameters  of  25  and  150  p  were  used  in  the 
configuration  perpendicular  to  the  vertical  (z)  direction  (see  Fig.  Bl). 

A  typical  curve  from  the  trough  measurements  is  given  in  Fig.  3,  both 
magnitude  and  phase.  The  frequency  along  the  x  axis  is  that  of  the  acoustic 
field.  The  recorded  frequency  is  twice  this  value.  This  should  be  kept  in 
mind  in  observing  e.g.  electrically  picked-up  signals  like  those  from  the  ac 
power  lines.  There  is  practically  no  phase  shift  between  current  through  the 
coil  and  output  from  the  pressure  hydrophone.  Thus  the  recorded  phase  from 
the  anemometer  is  directly  relative  to  the  particle  displacement.  The  pres¬ 
sure  field  is  approximately  independent  of  frequency  in  the  given  range.  Thus 
the  voltage  output  is  roughly  proportional  to  the  inverse  fourth  power  of  the 
frequency.  A  comparison  was  made  between  the  response  of  a  1211-20  sensor  and 
a  1210-20  sensor.  The  difference  is  in  the  orientation.  The  1211  sensor  has 
its  axis  vertical  (z  direction);  the  1210 's  is  horizontal  (y  direction). 

There  was  no  noticeable  difference  for  the  few  cases  where  these  were  compar¬ 
ed.  This  appears  to  indicate  that  there  is  no  strong  coupling  between  the 

heat  transfer  due  to  the  horizontal  motion  and  the  vertical  thermal  plume  (at 

-2 

least  at  the  given  Gr  -  6x10  ):  the  convection  is  different  for  a  vertical 

as  compared  with  a  horizontal  cylinder  (see  Ref.  26).  Shown  in  Table  3  is  the 
sensitivity  of  the  anemometer  to  rotation  of  the  sensor  in  the  plane  of  cylin¬ 
der  and  particle  motion  direction.  The  x  direction  is  the  direction  of  hori¬ 
zontal  particle  motion;  the  z  axis  is  vertical.  One  sees  that  there  is  a  very 
strong,  and  therefore  experimentally  useful,  directional  effect  at  low 
frequency.  It  decreases  drastically  with  increasing  frequency,  though. 


60 

Frequency  (Hi) 


The  frequency  of  the  abscissa  is  that  of  the  acoustic 
field.  The  received  frequency  is  twice  that  amount. 
The  upper  curve  shows  phase  and  the  curves  1,  2,  3  are 
magnitude  at  three  settings  of  the  analyzer: 

1  -  25  mV/division 

2  -  2.5  mV/division 

3  -  0.25  mV/division 


Sensor:  1211-20W;  0.  -  30°C;  3A  through  coil. 


Fig.  3  -  Example  of  anemometer  response  in  trough 


In  Table  4  the  velocities  are  given  that  refer  to  Tables  5  through  7. 

The  speeds  are  in  units  of  nun/s  and  are  independent  of  the  sensor  used,  since 
they  are  a  function  of  the  acoustic  field  only.  They  are  approximately 
proportional  to  the  inverse  of  the  frequency. 


Table  4  -  RMS  velocity  u(mm/s),  for  water  in  trough  as  a 
function  of  frequency  and  current  through 
shaker  coil.  (To  be  used  with  Tables  5-7.) 


Frequency 


(Hz) 

I  -  1A 

I  «  2A 

I  *  3A 

20 

2.60 

5.20 

7.80 

24 

2.30 

4.60 

6.90 

40 

1.37 

2.74 

4.71 

48 

1.27 

2.54 

3.82 

72 

0.713 

1.43 

2.14 

80 

0.656 

1.31 

1.97 

120 

0.449 

0.897 

1.35 

160 

0.372 

0.744 

1.12  1 

An  overview  of  the  data  for  water  is  given  in  Tables  5  through  7.  The 
dependent  variable  is  the  relative  displacement  C/d  computed  from  the  pressure 
hydrophone  calibration.  The  last  two  digits  10,  20,  and  60  in  the  sensor 
designation  indicate  diameters  of  25,  50,  and  150  pm,  respectively.  Three 
temperature  differences  Qg  between  sensor  and  medium  were  chosen:  15,  30,  and 
45°C.  The  upper  limit  for  usable  temperature  of  the  film  is  determined  by  the 
deleterious  influence  of  air  bubbles  driven  out  of  solution  and  adhering  to 
the  film  and  eventually  boiling.  The  dependence  of  the  output  ac  voltage  on 
the  relative  displacement  is  expected  to  be  quadratic,  and  b  is  the 
coefficient  in  the  quadratic  dependence  of  Nuac  (proportional  to  Vac)  on  the 
relative  displacement.  This  Is  further  discussed  in  Section  3.2  with  the  role 
of  the  Fourier  number,  which  is  an  entry  in  the  tables. 
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anemometer  output  from  horizontal  particle  m 
parentheses  ignored;  inconsistent. 
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Table  7  -  Hot-film  anemometer  output  from  horizontal  particl 
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The  variation  of  the  coefficient  b  with  frequency  and  sensor  prompted  the 
question  what  the  influence  of  the  viscosity  and  dlffusivlty  might  be. 
Therefore  the  experiments  were  repeated  with  ethylene  glycol  as  the  medium. 

It  has  a  thermal  diffuslvity  and  expansion  coefficient  comparable  to  those  of 
water,  but  the  viscosity  is  an  order  of  magnitude  larger.  The  ratio  of  the 
kinematic  viscosity  coefficients  at  40°C  is  13;  the  ratio  of  the  Prandtl 
numbers  at  40°C  is  21.  These  results  are  shown  in  Tables  9  through  11  while 
Table  8  gives  the  corresponding  speeds. 


Table  8  -  RMS  velocity  (mm/s)  for  ethylene 

glycol  in  trough  as  a  dimension  of 
frequency  and  current  through  shaker 
coil.  To  be  used  with  Tables  9-11. 


Hz 

I  =*  1A 

I  =  2A 

I  =  3A 

20 

2.58 

5.18 

7.70 

40 

1.47 

2.93 

4.40 

70 

0.776 

1.55 

2.32 

100 

0.578 

1.14 

1.73 

140 

0.453 

0.906 

1.36 

•.'•..'■i.VvV  v  -■ 
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anemometer  output  from  horizontal  particle  motion 


2. 12x10" J  87  140  0.00343  0.018  0.00687  0.025  0.0103  0.087  612  211 


horizontal  particle  motion 


Nu  *  1.33  100  70  0.0117  —  0.0235  0.162  0.0353 
Pr  -  95.3  101  100  0.00613  —  0.0123  0.037  0.0184 
Gf  -  8.41xl0-3  102  140  0.00343  —  0.00687  0.025  0.0103 
Nu,„/V_ _  -  0.466  V"1 


r  -V  \  1 


0. 401 


2.2.  Vertical  Particle  Motion 

The  vertical  particle  motion  is  created  in  an  NRL-USRD  type  G40 
calibrator.  The  response  of  the  anemometer  is  measured  at  the  same  frequency 
as  that  of  the  acoustic  field.  A  typical  curve  is  shown  in  Fig.  4.  Thanks  to 
the  greater  sensitivity  of  the  anemometer  in  this  configuration  it  was 
possible  to  obtain  useful  results  up  to  300  Hz.  Since  the  pressure  in  the 
calibrator  is  practically  independent  of  frequency,  the  response  is  roughly 
inversely  proportional  to  the  square  of  the  frequency.  The  difference  of  the 
anemometer  output  phase  and  the  phase  of  the  pressure  measurement  is 
practically  independent  of  frequency,  with  a  value  of  20®.  The  effect  of 
direction  on  sensitivity  was  measured  by  comparing  the  response  of  a  1210 
sensor  with  that  of  a  1211  sensor.  The  former  is  perpendicular  to  the 
particle  motion,  the  latter  is  parallel.  The  ratio  is  about  4  at  100,  200, 
and  300  Hz. 
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The  frequency  along  the  abscissa  and  that  of  the  anemometer  signal 
are  the  same,  equal  to  the  frequency  of  the  acoustic  field.  The 
upper  curve  shows  the  phase,  the  curves  1,  2,  3,  4  magnitude. 

1  -  25  mV  per  division 

2  -  10  mV  per  division 

3  -  2.5  mV  per  division 

4  -  1  mV  per  division 

Current  through  coil  is  400  mA.  Sensor  1210-20W,  ©s  *  45°C. 


Fig.  4.  Example  of  anemometer  response  in  G40  calibrator. 


CVaT* 


The  measurements  were  taken  for  water  only.  The  results  are  shown  in 
Tables  13  and  14.  The  two  sets  of  data  were  taken  about  a  year  apart. 

Table  12  gives  the  velocities  in  mm/s  corresponding  to  the  displacements  in 
Table  13.  The  dependence  of  the  ac  output  voltage  on  relative  displacement 
and  other  parameters  is  discussed  in  Section  3.3. 

Some  unusual  problems  were  encountered  during  the  experiments  in  the  G40 
calibrator.  Bubbles  were  formed  on  the  sensor,  resulting  in  irreproducible 
and  erratic  behavior.  This  same  type  of  problem  had  been  encountered  before. 
The  experiments  in  the  trough  had  not  shown  any  of  this  trouble.  It  was 
realized  that  the  only  difference  correlated  with  this  behavior  was  the  fact 
that  the  last  series  of  G40  experiments,  like  the  old  set,  was  done  after  the 
water  in  the  calibrator  had  been  exposed  to  wiping  off  the  walls  with  Aerosol 
0T100  surfactant  (American  Cyanamid  Corp.).  The  proper  operation  of  the 
anemometer  was  restored  only  after  cleaning  the  calibrator,  filling  with  fresh 
water,  spraying  the  sensors  with  methylalcohol ,  and  brushing  with  a  fine  paint 
brush.  A  full  explanation  was  not  attempted,  but  there  is  a  strong  suggestion 
that  the  surfactant,  while  aiding  the  release  of  a  bubble  that  has  been 
formed,  promotes  the  formation  of  small  bubbles  on  the  hot  sensor.  The  hot 
film  in  an  air-saturated  water  medium  creates  a  supersaturated  air  solution  by 
increasing  the  local  temperature.  Bubbles  may  not  be  formed  due  to  an  energy 
barrier  formed  by  the  minimum  bubble  size  that  is  energetically  possible  in 
connection  with  the  very  small  radius  of  the  sensor.  If  the  surfactant 
reduces  the  specific  surface  energy,  that  energy  barrier  may  be  lowered.  Even 
with  these  precautions  the  sensors  were  not  totally  reduced  to  their  old 
status.  Figure  5  shows  another  peculiar  effect  connected  with  this.  Before 
brushing  the  sensor,  the  curve  shows  a  very  regular  pattern  of  resonances 
above  about  100  Hz.  This  certainly  cannot  be  simple  mechanical  resonance  of 
bubbles  in  the  1/2-mm  size  and  smaller.  But  conceivably  it  is  an  effect  of 
one  or  a  whole  series  of  bubbles  that  are  compressed  and  expanded  by  the  film 
at  the  acoustic  frequency,  the  heat  transfer  changing  periodically.  Since 
heating  of  an  insulator  is  a  slow  process,  this  might  produce  a  low-frequency 
resonance.  Further  investigation  of  this  effect  would  be  of  interest. 


anemometer  output  from  vertical  particle  motion  in  water 

Data  of  Jul  &  Aug  1982 
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anemometer  output  from  vertical  particle  motion  in  water.  Data  of  Jul  1983 


(v  is 

vertical 

velocity) 


The  data  taken  in  the  G40  calibrator  when  presented  as  a  function  of  the 
Grashof  number  fall  on  two  different  lines  (see  Sec.  3.3).  This  might  be 
connected  with  the  same  bubble  problem.  No  surfactant  was  used  in  the  trough, 
and  no  bubble  trouble  was  encountered  there. 

2.3.  Imposed  Bias  Flow 

A  dc  bias  flow  was  imposed  by  means  of  a  nozzle  arrangement  attached  to  a 
yoke.  This  way  the  direction  of  the  water  jet  produced  by  the  nozzle  could  be 
swung  in  a  vertical  plane,  perpendicular  to  the  hot  film.  The  speed  of  the 
water  jet  was  regulated  by  a  needle  valve.  The  water  supply  to  the  nozzle  was 
gravity  driven  to  ensure  steady  laminar  flow.  The.  experiments  were  done  in 
the  G40  calibrator  where  the  particle  motion  is  vertical.  The  value  of  the 
bridge  voltage  of  the  anemometer  is  a  direct  indication  of  the  speed  of  the 
water  jet.  This  was  calibrated  in  a  rotating  tank.  The  relation  between  dc 
voltage  and  velocity  may  be  described  by  an  expression  of  the  type  given  in 
Eq.  (7)  or  (8).  The  diameter  of  the  nozzle  opening  was  3  mm,  and  its  front 
face  was  located  2  mm  from  the  hot  film.  The  results  have  a  semi-quantitative 
character  only;  more  data  will  be  taken  in  the  future.  The  conclusion  is 
warranted,  though,  that  the  sensitivity  increases  by  turning  on  the  jet.  With 
the  1210-20W  sensor,  0g  =  30°C,  and  a  jet  speed  of  3.3  cm/s,  an  increase  in 
voltage  output  is  observed  of  a  factor  3  at  frequencies  from  100  to  300  Hz  if 
the  jet  is  directed  vertically  upward.  When  the  jet  is  horizontal,  no  change 
in  ac  voltage  output  is  seen.  The  speed  of  3  cm/s  is  a  factor  10  above  the 
critical  speed  between  free  and  forced  convection,  from  Table  1.  Apparently 
the  free  convection  is  still  not  negligible,  since  directing  the  jet 
vertically  down  results  in  a  smaller  improvement  of  the  output  over  the  no- jet 
situation. 

3.  DATA  REDUCTION  AND  DISCUSSION  OF  RESULTS 

3.1.  Results  for  dc  Heat  Transfer 

It  is  of  interest  to  consider  the  relation  of  the  Nusselt  number  Nu^  to 
the  relevant  dimensionless  number(s).  Often  the  relation  is  given  [23]  in  a 
form  similar  to  Eq.  (9): 


B(GrPr)  . 


(13) 


Notice  that  the  product  GrPr  is  the  Rayleigh  number  and  may  be  considered  as  a 
Peclet  number  for  free  convection  if  Gr«l. 

A  linear  least-squares  analysis  was  performed  on  the  logarithmic  version 
of  Eq.  (13) 


log  =  logB  +  mlog(GrPr).  (14) 

The  results  of  this  analysis  for  the  three  sensors  show  that  the  1221-20 
sensor  does  not  fit  in  well  with  the  other  sensors.  This  can  be  understood  by 
noticing  that  the  1211  sensor  is  arranged  vertically  and  free  convection  from 
a  vertical  cylinder  is  different  from  that  of  a  horizontal  cylinder  [26]. 
Without  the  1211  data  one  finds  for  water  B  =  0.066  and  m  =  0.13  and  for 
ethylene  glycol  B  =  0.037  and  m  =  0.14. 

3.2.  Horizontal  Particle  Motion 

The  first  point  of  concern  is  to  check  the  expected  quadratic  dependence 
of  the  output  voltage  V ac  on  the  relative  displacement  5/d.  The  relative 
displacement  at  given  frequency  and  temperature  difference  is  proportional  to 
the  current  through  the  shaker  coil.  This  current  was  varied:  1,  2,  and  3A. 
All  through  Tables  5,  6,  and  7  the  quadratic  relation  is  quite  well  followed. 
There  is  some  decrease  at  higher  values  of  5/d,  at  about  2.5,  which  might  in¬ 
dicate  the  appearance  of  the  next  higher  term,  a  fourth  power  in  5/d.  The 
coefficient  b  in  the  relation  between  the  rms  Nusselt  number  Nuac  and 
displacement 


Nuac  =  b(5/d)2  (15) 

2 

was  found  by  averaging  V  /(5/d)  for  the  three  current  values  and  multiplying 

dO 

the  results  by  the  ratio  Nuac/Vac  (see  Appendix  A).  This  amounts  to  a  least- 
squares  adjustment  whereby  the  observations  are  weighted  inversely  proportion¬ 
al  to  the  square  of  their  magnitude.  If  this  weighting  is  not  applied,  the 
larger  values  of  the  observations  overshadow  the  small  ones.  These  values  for 


b  are  given  in  the  next  to  last  column  in  Tables  5,  6,  and  7.  It  is  most 

striking  that  in  contrast  to  ihe  results  for  vertical  particle  motion,  this 

coefficient  b  is  dependent  on  frequency.  In  some  cases  there  is  little 

variation  with  o>,  see  data  sets  88-92,  93-97,  103-107,  &  108-112.  It  also 

varies  with  sensor  diameter  but  relatively  little  with  t  mperature  difference 

0g.  The  dimensionless  number  next  in  line  in  the  temperature  Eq.  (6)  is  the 

2 

Fourier  number  Fo  =  x/uxi  .  The  values  of  the  inverse  Fourier  number  are  given 
in  the  last  column  to  the  right.  There  is  an  obvious  correlation  between  the 
variation  of  b  and  the  increase  of  Fo  To  study  this  quantitatively  a 
least-squares  regression  was  performed  on  a  relation  between  b  and  the  Fourier 
number  of  the  form 


log  b  =  log  c  -  mlog  Fo. 


(16) 


The  results  are  shown  in  Table  15  for  each  sensor  at  given  0  .  From  the 

s 

least-squares  regression  one  may  derive  an  rms  residual  r  according  to  the 
expression 


2  .  2  (1°8bm  "  log  bc)2 
N  -  2 


(17) 


where  N  is  the  number  of  observations  in  the  regression,  bm  is  the  measured 

value,  and  b  is  the  calculated  value.  An  "rms  ratio"  R.  „  may  be  connected 

with  r  by  log  Rrms  =  t.  This  ratio  is  used  for  comparison  of  the  various 

regression  models  as  a  test  on  the  goodness  of  fit.  The  results  show  that  the 

exponent  m  is  reasonably  constant,  but  the  coefficient  c  varies  with  the 

sensor  diameter.  Thus,  the  scaling  with  w  and  d  is  not  compatible  with 
2 

a  d  w  combination  as  implied  by  Fo  or  St  scaling.  If  one  pools  the  values  for 
the  same  sensor  at  different  Qg ,  one  finds  the  following  results.  For  a  water 

medium:  c  »  9.8x10  ,  m  »  0.55,  and  Rrms  *  1.2  for  the  25-pm  sensor; 

-3  -3 

c  ■  16.6x10  ,  m  -  0.69,  and  R^g  =  1.2  for  the  50-pm  sensor;  and  c  *  64xio  , 

m  ■  0.46,  and  Rrms  ”  1.4  for  the  150-pm  sensor.  In  the  case  of  glycol  the 

_3 

corresponding  numbers  are:  c  =  5.9x10  ,  m  »  0.31,  and  Rrms  ”  1.4  for  the 

25-um  sensor;  c  *  11.8x10  ,  m  «  0.34,  and  Rrms  *  1.4  for  the  50-pm  sensor; 

-  7.4xlO-3,  ra  - 


and  c 


0.67,  and  R 


1.6  for  the  150-pm  sensor.  All  the 


water  data  together  give  c  =  lO.AxlO-^,  m  =  0.90,  and  Rrms  =  1.4.  All  the 
glycol  data  together  give  c  =  4.5x10'^,  m  =  0.75,  and  Rrms  =  1.6. 

The  results  for  the  coefficient  c  and  exponent  m  for  ethylene  glycol 
differ  considerably  from  those  for  water.  This  suggests  that  a  factor  other 
than  the  Fourier  number  is  involved.  Since  the  viscosity  of  ethylene  glycol 
is  much  larger  than  that  of  water,  one  might  assume  that  a  combination  of  Pr 
and  Gr  would  provide  a  better  fit,  in  cooperation  with  the  Fo  regression. 
Therefore,  a  least-squares  adjustment  was  performed  on  the  expression 


log  b  =  log  c  +  m^logx^  +  m^logx^ ,  (18) 

where  the  first  variable  x^  may  be  chosen  to  be  either  Fo-^  or  St~^,  while  the 

second  variable  X2  may  be  equated  with  Pr,  Gr,  or  Ra(=GrPr).  The  involvement 

of  Gr  promises  to  rectify  in  part  the  scaling  with  respect  to  d  and  00,  which 

2 

was  poorly  represented  by  the  combination  d  u>.  Out  of  the  six  possible 
combinations  a  definite  preference  was  shown  for  Fo  and  Gr,  judging  by  the  rms 
ratio,  by  the  randomness  of  the  measured  value  bm  divided  by  the  computed 
value  bc,  and  by  the  standard  deviations  and  cross  correlations  between  the  c, 
m^ ,  and  m2*  Thus,  the  best  fit  appears  to  be  given  by  the  following 
expression  for  b: 


b  =  0.042  Fo~°’ 56  Gr0'22. 


The  results  for  all  the  values  of  b  are  given  in  Table  16.  Of  course,  the 

0.22 

product  0.042  Gr  ‘  is  equivalent  to  the  former  (variable)  coefficient  c  of 
the  regression  on  Fo.  This  product  is  listed  for  comparison  in  the  last 
column  of  Table  15. 

The  values  and  standard  deviations  are  (1)  logc  =  -1.38  ±0.03, 

(2)  m^  =  0.56  ±0.02,  (3)  m2  =  0.22  ±0.01  and  the  cross  correlations  are 

P^2  =  -0.80,  =  0.83,  P23  *  -0.50.  The  scaling  in  terms  of  d  and  u>  implied 

by  Eq.  (19)  is  approximately  d2u>  ^  . 

The  dependence  of  the  coefficient  b  on  the  Fourier  and  Grashof  numbers 
means  that,  given  a  constant  ratio  5/d,  the  sensor  is  more  efficient  for 
larger  Grashof  number  and  smaller  Fourier  number.  Physically  one  might 
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regression  of  b  with  respect 
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tentatively  interpret  this  as  follows.  The  inverse  Fourier  number  expresses 
the  ratio  of  local  heating  to  heat  transport  by  conduction.  If  this  ratio  is 
small,  the  acoustically  induced  heating  and  cooling  are  obviously  small 
compared  with  the  average  heat  conduction.  The  dependence  on  the  Grashof 
number  suggests  that  for  small  values  of  Gr  the  hot  film  is  surrounded  by  a 
viscosity  dominated  region  extending  many  diameters  from  the  sensor,  on  which 
the  variable  acoustic  flow  has  little  influence.  For  larger  Gr  the  viscous 
region  is  reduced  to  a  boundary-type  layer,  and  the  cooling  effect  of  the 
acoustic  flow  can  approach  the  hot  film  more  closely. 

The  pattern  of  phase  variation  with  frequency  (Fig.  3)  also  changes  with 
d  and  oj ;  a  complete  study  of  its  scaling  properties  has  yet  to  be  done.  It  is 
possible  that  a  simpler  scaling  law  could  be  discovered  for  in-phase  and  out- 
of-phase  components  separately  (with  respect  to  the  phase  of  the  pressure  or 
displacement  in  the  main  body  of  the  medium).  Of  course,  this  could  not  be 
used  for  measurement  of  an  unknown  field  since  the  absolute  phase  would  not  be 
known,  but  it  might  be  useful  for  better  insight  into  the  heat  transfer 
process. 


Calculated  value  b.  and  ratio  of  measured  to  calculated  value 
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3.3.  Vertical  Particle  Motion 


Inspection  of  Tables  13  and  14  shows  that  Nuac  is  quite  well  represented 
by  a  linear  relationship  with  respect  to  the  relative  displacement,  since  it 
varies  proportionally  to  the  current  through  the  calibrator  coil.  The 
coefficient  a  in  the  relation 

Nuac  -  a(c/d)  (20) 

was  computed  by  averaging  Nu  /(5/d)  for  the  various  current  values.  This 

ac 

amounts  to  a  least-squares  regression  whereby  the  observed  values  are  weighted 
by  a  factor  inversely  proportional  to  the  magnitude  squared  of  the 
measurement. 

It  is  seen  that  the  value  of  a  is  reasonably  independent  of  frequency 
within  the  accuracy  limits  of  the  experiment.  Since  the  ac  particle  motion  is 
superimposed  on  the  free  convection,  it  is  logical  to  attempt  a  regression  of 
the  coefficient  a  on  the  Grashof  or  Rayleigh  number.  Figure  6  shows  the 
relationship  for  the  two  sets  of  data  in  1982  and  1983.  It  Is  quite  striking 
that  the  data  appear  to  fall  on  two  distinct  lines  not  related  to  sensor  type 
or  time  of  measurement.  No  explanation  for  this  unfortunate  discrepancy  has 
been  found  to  date.  It  is  possible  that  the  formation  of  bubbles,  as 
described  in  Section  2.2,  played  a  role;  but  it  is  difficult  to  see  how  this 
would  have  created  two  distinct  lines  and  not  a  continuum,  of  values.  A  least- 
squares  regression  was  performed  on  the  data  of  the  two  separate  cases  with 
the  logarithmic  dependence 

log  a  -  log  c  +  m^log  Ra.  (21) 

The  fit  for  the  two  sets  was  quite  good.  The  values  for  the  coefficient  cr^ 
and  the  exponent  m^  were  0.217  and  0.360  with  R^g  “  1.04  for  the  upper  line 
and  0.130  and  0.315  with  Rrms  *  1.06  for  the  lower  line.  The  data  coverage 
was  insufficient  to  distinguish  between  Gr  and  Ra  as  the  better  independent 
variable.  To  resolve  this  question,  measurements  with  ethylene  glycol  in  the 
calibrator  would  be  in  order. 


The  upper  number  between  parentheses  is  the  temperature  difference  0 
The  lower  number  is  the  sensor  diameter  in  p. 

•  -  data  of  Jul  83 
■  -  data  of  Jul-Aug  82. 


6  -  Coefficient  of  relation  Nuflc 
number. 


a(C/d)  as  a  function  of  Grashof 


4.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 


The  present  study  shows  that  hot-film  anemometry  can  be  successfully 
applied  to  the  detection  of  particle  motion  in  hydroacoustic  fields.  Its 
sensitivity  at  constant  intensity  increases  for  decreasing  frequency,  in 
contrast  to  velocity  hydrophones  based  on  accelerometers.  It  is  felt  that, 
conversely,  the  anemometer  response  to  harmonic  acoustic  fields  may  also  be 
used  as  a  tool  to  study  free  convection  about  a  cylinder. 

The  prime  area  to  be  developed  is  the  detection  of  particle  motion  under 
imposed  bias  flow.  It  has  greater  sensitivity,  frees  the  experiment  of  the 
dependence  on  the  direction  of  gravity,  and  adds  a  valuable  control  on 
determination  of  direction  of  the  acoustic  motion. 

The  discrepancy  in  the  data  in  vertical  particle  motion  should  be 
investigated. 

Further  investigation  of  the  empirical  relations  in  terms  of 
dimensionless  groups  is  in  order,  including  also  the  variation  in  phase. 

Overall  improvement  of  the  accuracy  of  the  measurements,  which  at  this 
time  is  estimated  at  about  10  to  20%,  appears  usaful.  This  would  include 
modifications  of  trough  and  vertical  calibrator. 
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Appendix  A 


COMPUTATION  OF  DIMENSIONLESS  NUMBERS 

Most  of  the  data  on  the  properties  of  glycol  and  water  were  taken  from 
Ref.  Al.  The  data  for  the  expansion  coefficient  of  ethylene  glycol  are  found 
in  Ref.  A2.  The  data  near  40°C  were  interpolated  by  a  second-degree 
polynomial  through  the  data  points  at  20,  40,  and  60#C.  T  is  in  degree 
Celsius. 

For  water:  k(W  m_1K_1)  =  0.552  +  2.35xlO~3T  -  11.2xl0~6T2 

v(m2s“l)  =  1.522xlO-6  -  30.0xl0_9T  +  0.210xlO_9T2 
<x(K-1)  =  -12.87xl0-6  +  11.99x10-6T  -  0.0509xl0_6T2 
Pr  =  11.06  -  0.236T  +  1.7xlO_3T? 

For  ethylene  glycol: 

k(Wm_1K_1)  =  0.239  +  0.575xl0~3T  -  3.75xlO_6T2 
v(m2s_1)  =  36.23xlO~6  -  1.016xl0~6T  +  8.19xl0_9T2 
a(K_1)  -  0. 556*10“ 3  +  8.70xl0_6T  -  0.145xl0_6T2 
Pr  *  384  -  10. 7T  +  86.3xl0~3T? 

The  material  parameters  are  computed  at  an  average  temperature  of  the 
sensor  and  medium  temperatures  (Tg  +  TQ)/2.  The  Nusselt  number  is  defined  as 
Hd/kQ8,  where  H  is  the  heat  flux.  H  is  computed  from  the  measured  bridge 
voltage  in  the  following  way.  The  sensor  resistance  R  is  in  series  with  the 
bridge  resistor  R^  in  the  anemometer  bridge,  with  a  value  Rb  =  40  fl.  The 
resistance  measured  in  the  bridge  R^  includes  the  resistance  of  the  cable, 
probe-holder,  and  probe  in  addition  to  the  sensor  resistance  itself.  The 
cable  resistance  is  0.1  ft,  that  of  the  probe-holder  is  0.08  ft,  and  the 
resistance  Rp  of  the  probe  varies  with  the  given  sensor. 

• 

The  heat  transfer  rate  Q  is  given  by 


Q 


2 

dc 


(R,+V 


2* 


(Al) 


For  a  typical  measurement  with  a  1210-20W  sensor  with  0g  =  30°,  R  *  6  (!, 

one  finds  V.  =  4.78  V.  Thus  5  =  64  mW.  The  heat  flux  is  determined  by 
dc  x 

H  =  Q/n£d  (A2 

and  the  dc  Nusselt  number  Nucjc  by 

““dc  ■  ''dcR/’*kVRt  +  <« 

The  factor  needed  to  convert  amplitude  or  rms  value  of  the  ac  voltage  to  the 

corresponding  value  for  the  Nusselt  number  variation  is  given  by 

Nu  /V  =  2Nu ,  /V_,  . 
ac  ac  dc  dc 
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TROUGH  AND  SHAKER  ARRANGEMENT 


The  horizontal  particle  movement  is  produced  by  a  trough,  especially 
designed  for  this  application,  attached  to  a  horizontal  shaker  (Fig.  Bl).  The 
original  trough  was  made  of  a  polyvinyl  chloride  (PVC)  cylinder  with  6.35-mm 
(1/4-in.)  wall  thickness  and  12.7-mm  ( 1/2-in. )-thick  lucite  endplates.  The 
PVC  trough  was  suspended  by  1. 5-m  (60-in.)  cables  from  the  ceiling.  An 
improved  version  was  constructed  of  steel.  Its  dimensions  are  shown  in  cross 
section  in  Fig.  B2.  Notice  that  the  trough  has  a  free  surface  of  the 
liquid.  A  similar  device  was  described  by  Bauer  [Bl],  based  on  a  closed 
cylinder.  A  later  version  has  a  free  surface  [ B2 ] .  The  steel  trough  was 
suspended  from  an  A-frame  with  four  steel  wires  of  0.76-mm  (0.030-in.)  diam. 
and  20-cm  length. 


Fig.  Bl  -  Trough  for  horizontal  particle  motion 


Analysis  of  the  resulting  motion  is  more  complicated  as  a  consequence  of 
the  free  surface,  especially  for  a  cylindrical  geometry.  One  may  readily 
solve  the  equation  of  motion  under  the  proper  boundary  conditions  for  a 
rectangular  trough.  Although  the  frequencies  are  rather  low  in  the  present 
study,  it  is  assumed  for  generality  that  the  pressure  is  a  solution  of 
Helmholtz's  equation 

(V2  +  ko2)p  =  0,  (Bl) 

where  kQ  is  the  wavenumber  in  the  medium  kQ  =  co/c0 ,  with  cQ  the  propagation 
speed.  As  to  the  boundary  conditions,  assuming  rigidity  of  the  trough,  the 
end  faces  at  x  =  ±L/2  impose  a  condition  on  the  velocity 

u(-L/2  ,  z)  =  u(L/2,z)  =  U  (B2) 

and  the  vertical  velocity  w  is  zero  at  the  bottom  and  thus  the  pressure 
gradient 


^  *  0  at  z  =  0.  (B3) 

It  is  assumed  that  the  surface  variations  are  not  important  dynamically,  thus 
diverging  from  the  usual  development  of  surface  waves,  to  give 

p(x,h)  =  0,  (B4) 


where  h  is  the  water  depth.  The  solution  for  pressure  and  velocity  components 
is 
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The  trough  Is  connected  to  a  horizontal  shaker.  The  pressure  is  measured 
by  a  USRD  type  F61  standard  hydrophone  [B3].  Figure  B3  shows  the  pressure 
variation  along  the  axis  of  the  PVC  trough  starting  from  the  side  attached  to 
the  shaker.  The  solid  curve  is  a  fourth-degree  polynomial  fit. 

Figure  B4  gives  the  pressure  profile  for  the  steel  trough.  The  symmetry 
of  the  curve  about  the  center  of  the  trough  has  obviously  improved  thanks  to 
the  greater  rigidity  of  the  steel.  Moreover,  the  slit  in  the  top  of  the 
cylinder  was  smaller  than  for  the  PVC  trough.  Thus  the  free  surface  was 
smaller  than  in  that  case.  As  a  consequence  the  vertical  pressure  gradient 
needed  to  accelerate  fluid  towards  the  free  surface  is  less.  Therefore  the  S 


is  less  deep  and  more  like  the  straight  line  that  would  prevail  for  a 
completely  closed  cylinder.  The  polynomial  fit  served  to  find  the  derivative 
in  the  center  of  the  trough  where  the  hot-film  anemometer  measurements  were 
performed.  A  set  of  derivative  curves  at  various  frequencies  is  shown  in 
Fig.  B5  in  terms  of  the  relative  displacement  (5/d). 


V  V.-  .H  .  .  .T  . 


DISTANCE  (cm) 


Fig.  B5  -  Relative  displacement  in  steel  trough  as  a  function  of 
distance  along  axis  for  various  frequencies;  d  ■  50p 
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